1 BCO comparison of CASAL and Casal2 model configurations

This document compares the results of at least 2 CASAL model configurations (base and at least one sensitivity) and up to 8 Casal2 model configurations (3 BetaDiff, 2 CppAD, and 3 ADOL-C).

The CASAL model sensitivity 1 has a smaller minimisation tolerance value than the CASAL base model (1e-7 vs. 2e-5).

The Casal2 ADOL-C and BetaDiff low tolerance models have a smaller tolerance value than the CASAL base model (1e-7 vs. 2e-5). The Casal2 CppAD models have a tolerance value of 1e-9.

1.1 BCO model characteristics

The main characteristics of the Test Case BCO (ling) CASAL model are:

  • one stock, ages 3 - 20
  • one area
  • years 1900 - 2019, projection year 2020
  • one time step
  • two stock categories: male and female
  • natural mortality-at-age is 0.14 for males and 0.20 for females
  • sex-specific logistic maturity
  • sex-specific von Bertalanffy age-length relationship
  • one length-weight relationship (\(W = aL^b\))
  • Beverton-Holt stock-recruitment relationship, with steepness (h) 0.75 and \(\sigma_R\) 0.6
  • no ageing error
  • one survey with sex-specific age-based logistic selectivity
  • three fisheries, all with sex-specific length-based logistic selectivity and retained/discarded catches

Observation data include:

  • survey index
  • fishery CPUE index
  • fishery proportions-at-length

Parameters estimated include:

  • B0
  • male \(L_{\infty}\)
  • male age-length CV2
  • sex-specific survey logistic selectivity (a50 and ato95)
  • YCS (recruitment deviations)

1.2 R environment

## [1] "Mon Nov 23 15:38:54 2020"
## R version 3.6.0 (2019-04-26)
## Platform: x86_64-redhat-linux-gnu (64-bit)
## Running under: CentOS Linux 7 (Core)
## 
## Matrix products: default
## BLAS/LAPACK: /usr/lib64/R/lib/libRblas.so
## 
## locale:
##  [1] LC_CTYPE=en_NZ.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_NZ.UTF-8        LC_COLLATE=en_NZ.UTF-8    
##  [5] LC_MONETARY=en_NZ.UTF-8    LC_MESSAGES=en_NZ.UTF-8   
##  [7] LC_PAPER=en_GB.UTF_8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_NZ.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] casal2_1.0     casal_2.30     ggthemes_4.2.0 ggmcmc_1.5.0   tidyr_1.1.2   
##  [6] ggplot2_3.3.2  gridExtra_2.3  coda_0.19-4    rlist_0.4.6.1  huxtable_5.1.1
## [11] dplyr_1.0.2    plyr_1.8.6     devtools_2.3.2 usethis_1.6.3 
## 
## loaded via a namespace (and not attached):
##  [1] tidyselect_1.1.0   xfun_0.19          remotes_2.2.0      purrr_0.3.4       
##  [5] lattice_0.20-38    colorspace_1.4-1   vctrs_0.3.5        generics_0.1.0    
##  [9] testthat_2.3.2     htmltools_0.5.0    yaml_2.2.1         rlang_0.4.8       
## [13] pkgbuild_1.1.0     pillar_1.4.7       glue_1.4.2         withr_2.3.0       
## [17] RColorBrewer_1.1-2 sessioninfo_1.1.1  lifecycle_0.2.0    stringr_1.4.0     
## [21] munsell_0.5.0      gtable_0.3.0       memoise_1.1.0      evaluate_0.14     
## [25] GGally_1.4.0       knitr_1.30         callr_3.5.1        ps_1.4.0          
## [29] fansi_0.4.1        Rcpp_1.0.5         scales_1.1.1       desc_1.2.0        
## [33] pkgload_1.1.0      fs_1.5.0           digest_0.6.27      stringi_1.5.3     
## [37] processx_3.4.4     rprojroot_2.0.2    grid_3.6.0         cli_2.2.0         
## [41] tools_3.6.0        magrittr_2.0.1     tibble_3.0.4       crayon_1.3.4      
## [45] pkgconfig_2.0.3    ellipsis_0.3.1     data.table_1.13.2  prettyunits_1.1.1 
## [49] reshape_0.8.8      assertthat_0.2.1   rmarkdown_2.5      R6_2.5.0          
## [53] compiler_3.6.0

1.3 CASAL and Casal2 model output

source('../../R-functions/report_read_in_CASAL_MPD_files.R')
source('../../R-functions/report_read_in_Casal2_MPD_files.R')
## Warning in sqrt(1/diag(V)): NaNs produced
## Warning in cov2cor(data): diag(.) had 0 or NA entries; non-finite result is
## doubtful
## Warning in sqrt(1/diag(V)): NaNs produced
## Warning in cov2cor(data): diag(.) had 0 or NA entries; non-finite result is
## doubtful

1.4 Tables

Tables of parameter estimates and objective function components for the CASAL and Casal2 model MPD results

CASAL parameter estimates
rownamesBase_ModelSensitivity_1Percent_Diff
initialization.B0       14170.0000014170.000000.000
recruitment.YCS1     1.000001.000000.000
recruitment.YCS2     1.000001.000000.000
recruitment.YCS3     1.000001.000000.000
recruitment.YCS4     1.000001.000000.000
recruitment.YCS5     1.000001.000000.000
recruitment.YCS6     1.000001.000000.000
recruitment.YCS7     1.000001.000000.000
recruitment.YCS8     1.000001.000000.000
recruitment.YCS9     1.000001.000000.000
recruitment.YCS10    1.000001.000000.000
recruitment.YCS11    1.000001.000000.000
recruitment.YCS12    1.000001.000000.000
recruitment.YCS13    1.000001.000000.000
recruitment.YCS14    1.000001.000000.000
recruitment.YCS15    1.000001.000000.000
recruitment.YCS16    1.000001.000000.000
recruitment.YCS17    1.000001.000000.000
recruitment.YCS18    1.000001.000000.000
recruitment.YCS19    1.000001.000000.000
recruitment.YCS20    1.000001.000000.000
recruitment.YCS21    1.000001.000000.000
recruitment.YCS22    1.000001.000000.000
recruitment.YCS23    1.000001.000000.000
recruitment.YCS24    1.000001.000000.000
recruitment.YCS25    1.000001.000000.000
recruitment.YCS26    1.000001.000000.000
recruitment.YCS27    1.000001.000000.000
recruitment.YCS28    1.000001.000000.000
recruitment.YCS29    1.000001.000000.000
recruitment.YCS30    1.000001.000000.000
recruitment.YCS31    1.000001.000000.000
recruitment.YCS32    1.000001.000000.000
recruitment.YCS33    1.000001.000000.000
recruitment.YCS34    1.000001.000000.000
recruitment.YCS35    1.000001.000000.000
recruitment.YCS36    1.000001.000000.000
recruitment.YCS37    1.000001.000000.000
recruitment.YCS38    1.000001.000000.000
recruitment.YCS39    1.000001.000000.000
recruitment.YCS40    1.000001.000000.000
recruitment.YCS41    1.000001.000000.000
recruitment.YCS42    1.000001.000000.000
recruitment.YCS43    1.000001.000000.000
recruitment.YCS44    1.000001.000000.000
recruitment.YCS45    1.000001.000000.000
recruitment.YCS46    1.000001.000000.000
recruitment.YCS47    1.000001.000000.000
recruitment.YCS48    1.000001.000000.000
recruitment.YCS49    1.000001.000000.000
recruitment.YCS50    1.000001.000000.000
recruitment.YCS51    1.000001.000000.000
recruitment.YCS52    1.000001.000000.000
recruitment.YCS53    1.000001.000000.000
recruitment.YCS54    1.000001.000000.000
recruitment.YCS55    1.000001.000000.000
recruitment.YCS56    1.000001.000000.000
recruitment.YCS57    1.000001.000000.000
recruitment.YCS58    1.000001.000000.000
recruitment.YCS59    1.000001.000000.000
recruitment.YCS60    1.000001.000000.000
recruitment.YCS61    1.000001.000000.000
recruitment.YCS62    1.000001.000000.000
recruitment.YCS63    1.000001.000000.000
recruitment.YCS64    1.000001.000000.000
recruitment.YCS65    1.000001.000000.000
recruitment.YCS66    1.000001.000000.000
recruitment.YCS67    1.000001.000000.000
recruitment.YCS68    1.000001.000000.000
recruitment.YCS69    1.000001.000000.000
recruitment.YCS70    1.000001.000000.000
recruitment.YCS71    1.000001.000000.000
recruitment.YCS72    1.000001.000000.000
recruitment.YCS73    1.000001.000000.000
recruitment.YCS74    1.000001.000000.000
recruitment.YCS75    1.000001.000000.000
recruitment.YCS76    1.000001.000000.000
recruitment.YCS77    1.000001.000000.000
recruitment.YCS78    1.000001.000000.000
recruitment.YCS79    1.000001.000000.000
recruitment.YCS80    1.000001.000000.000
recruitment.YCS81    1.000001.000000.000
recruitment.YCS82    1.000001.000000.000
recruitment.YCS83    1.000001.000000.000
recruitment.YCS84    0.234800.234800.000
recruitment.YCS85    0.231600.231600.000
recruitment.YCS86    0.232200.232200.000
recruitment.YCS87    0.241200.241200.000
recruitment.YCS88    0.265100.265100.000
recruitment.YCS89    0.317200.317200.000
recruitment.YCS90    0.422400.422400.000
recruitment.YCS91    0.613900.613900.000
recruitment.YCS92    0.923000.923000.000
recruitment.YCS93    1.429001.429000.000
recruitment.YCS94    1.480001.480000.000
recruitment.YCS95    1.009001.009000.000
recruitment.YCS96    0.846200.846200.000
recruitment.YCS97    0.855500.855500.000
recruitment.YCS98    1.087001.087000.000
recruitment.YCS99    1.530001.530000.000
recruitment.YCS100   1.027001.027000.000
recruitment.YCS101   0.799900.799900.000
recruitment.YCS102   0.820800.820800.000
recruitment.YCS103   0.631800.631800.000
recruitment.YCS104   0.805200.805200.000
recruitment.YCS105   0.745300.745300.000
recruitment.YCS106   0.463400.463400.000
recruitment.YCS107   1.075001.075000.000
recruitment.YCS108   0.899200.899200.000
recruitment.YCS109   0.385100.385100.000
recruitment.YCS110   0.912700.912700.000
recruitment.YCS111   0.709200.709200.000
recruitment.YCS112   0.829000.829000.000
recruitment.YCS113   0.423300.423300.000
recruitment.YCS114   1.278001.278000.000
recruitment.YCS115   0.688100.688100.000
recruitment.YCS116   0.742500.742500.000
recruitment.YCS117   0.570500.570500.000
recruitment.YCS118   0.451600.451600.000
recruitment.YCS119   1.000001.000000.000
recruitment.YCS120   1.000001.000000.000
size_at_age.Linf_male47.1200047.120000.000
size_at_age.cv2_male 0.062560.062560.000
selectivity[potSsel].male1    3.039003.039000.000
selectivity[potSsel].male2    0.190700.190600.052
selectivity[potSsel].female1  3.316003.316000.000
selectivity[potSsel].female2  1.045001.045000.000
Casal2 parameter estimates
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_tolcppad_casal_flags_oncppad_casal_flags_offadolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
process[Recruitment].b0              14199.6000014210.7000014199.6000014199.6000014210.7000014199.6000014210.7000014199.60000
age_length[asMm0].linf            47.3306047.3135047.3306047.3306047.3135047.3306047.3135047.33060
age_length[asMm0].cv_last         0.059780.061420.059780.059780.061420.059780.061420.05978
selectivity[potSurveySel_male].a50             3.039353.040113.039323.039593.040153.039353.040113.03932
selectivity[potSurveySel_male].ato95           0.190780.191760.190620.191910.191970.190780.191760.19062
selectivity[potSurveySel_female].a50             3.314493.315093.314503.314503.315103.314493.315093.31450
selectivity[potSurveySel_female].ato95           1.033241.026471.033251.033251.026511.033241.026471.03325
process[Recruitment].ycs_values{1980}0.235220.235120.235220.235220.235120.235220.235120.23522
process[Recruitment].ycs_values{1981}0.232180.232090.232180.232180.232090.232180.232090.23218
process[Recruitment].ycs_values{1982}0.233210.233170.233210.233210.233170.233210.233170.23321
process[Recruitment].ycs_values{1983}0.241990.241990.241990.241990.241990.241990.241990.24199
process[Recruitment].ycs_values{1984}0.264640.264700.264640.264640.264700.264640.264700.26464
process[Recruitment].ycs_values{1985}0.314840.314950.314840.314840.314950.314840.314950.31484
process[Recruitment].ycs_values{1986}0.416170.416350.416170.416170.416350.416170.416350.41617
process[Recruitment].ycs_values{1987}0.601740.601990.601740.601740.602000.601740.601990.60174
process[Recruitment].ycs_values{1988}0.902120.902570.902120.902120.902550.902120.902570.90212
process[Recruitment].ycs_values{1989}1.378331.379871.378331.378331.379861.378331.379871.37833
process[Recruitment].ycs_values{1990}1.560021.558561.560011.560011.558581.560021.558561.56001
process[Recruitment].ycs_values{1991}1.028901.028501.028911.028911.028511.028901.028501.02891
process[Recruitment].ycs_values{1992}0.847880.847820.847880.847880.847800.847880.847820.84788
process[Recruitment].ycs_values{1993}0.846860.847080.846850.846850.847080.846860.847080.84685
process[Recruitment].ycs_values{1994}1.065361.065751.065361.065361.065751.065361.065751.06536
process[Recruitment].ycs_values{1995}1.548971.548051.548971.548971.548101.548971.548051.54897
process[Recruitment].ycs_values{1996}1.053291.053451.053271.053271.053371.053291.053451.05327
process[Recruitment].ycs_values{1997}0.803050.802900.803060.803060.802940.803050.802900.80306
process[Recruitment].ycs_values{1998}0.816690.816580.816690.816690.816580.816690.816580.81669
process[Recruitment].ycs_values{1999}0.628380.628940.628380.628380.628940.628380.628940.62838
process[Recruitment].ycs_values{2000}0.802060.802340.802060.802060.802340.802060.802340.80206
process[Recruitment].ycs_values{2001}0.745810.745990.745810.745810.745990.745810.745990.74581
process[Recruitment].ycs_values{2002}0.467130.466750.467130.467130.466750.467130.466750.46713
process[Recruitment].ycs_values{2003}1.094231.090631.094231.094231.090631.094231.090631.09423
process[Recruitment].ycs_values{2004}0.921280.917110.921280.921280.917100.921280.917110.92128
process[Recruitment].ycs_values{2005}0.384130.382560.384130.384130.382560.384130.382560.38413
process[Recruitment].ycs_values{2006}0.905890.901900.905890.905890.901900.905890.901900.90589
process[Recruitment].ycs_values{2007}0.709810.720450.709810.709810.720450.709810.720450.70981
process[Recruitment].ycs_values{2008}0.835570.834790.835570.835570.834790.835570.834790.83557
process[Recruitment].ycs_values{2009}0.426300.426060.426290.426300.426060.426300.426060.42629
process[Recruitment].ycs_values{2010}1.287461.287071.287461.287461.287071.287461.287071.28746
process[Recruitment].ycs_values{2011}0.690330.691100.690330.690330.691100.690330.691100.69033
process[Recruitment].ycs_values{2012}0.736380.736780.736380.736380.736780.736380.736780.73638
process[Recruitment].ycs_values{2013}0.561560.561900.561560.561560.561900.561560.561900.56156
process[Recruitment].ycs_values{2014}0.439780.439780.439780.439780.439780.439780.439780.43978
Casal2 parameter estimates: Percent Difference from betadiff_casal_flags_on
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_tolcppad_casal_flags_oncppad_casal_flags_offadolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
process[Recruitment].b0              0.000-0.0780.0000.000-0.0780.000-0.0780.000
age_length[asMm0].linf            0.0000.0360.0000.0000.0360.0000.0360.000
age_length[asMm0].cv_last         0.000-2.7570.0000.000-2.7570.000-2.7570.000
selectivity[potSurveySel_male].a50             0.000-0.0250.001-0.008-0.0260.000-0.0250.001
selectivity[potSurveySel_male].ato95           0.000-0.5180.084-0.594-0.6260.000-0.5180.084
selectivity[potSurveySel_female].a50             0.000-0.018-0.000-0.000-0.0180.000-0.018-0.000
selectivity[potSurveySel_female].ato95           0.0000.655-0.001-0.0010.6510.0000.655-0.001
process[Recruitment].ycs_values{1980}0.0000.0450.0000.0000.0450.0000.0450.000
process[Recruitment].ycs_values{1981}0.0000.0370.0000.0000.0370.0000.0370.000
process[Recruitment].ycs_values{1982}0.0000.0190.0000.0000.0190.0000.0190.000
process[Recruitment].ycs_values{1983}0.000-0.0020.000-0.000-0.0020.000-0.0020.000
process[Recruitment].ycs_values{1984}0.000-0.0210.0000.000-0.0210.000-0.0210.000
process[Recruitment].ycs_values{1985}0.000-0.0350.0000.000-0.0360.000-0.0350.000
process[Recruitment].ycs_values{1986}0.000-0.0420.0000.000-0.0420.000-0.0420.000
process[Recruitment].ycs_values{1987}0.000-0.042-0.000-0.000-0.0430.000-0.042-0.000
process[Recruitment].ycs_values{1988}0.000-0.0490.0010.001-0.0470.000-0.0490.001
process[Recruitment].ycs_values{1989}0.000-0.1120.0000.000-0.1110.000-0.1120.000
process[Recruitment].ycs_values{1990}0.0000.0940.0010.0010.0920.0000.0940.001
process[Recruitment].ycs_values{1991}0.0000.039-0.001-0.0010.0380.0000.039-0.001
process[Recruitment].ycs_values{1992}0.0000.0070.0000.0000.0090.0000.0070.000
process[Recruitment].ycs_values{1993}0.000-0.0270.0010.001-0.0260.000-0.0270.001
process[Recruitment].ycs_values{1994}0.000-0.0370.0000.000-0.0370.000-0.0370.000
process[Recruitment].ycs_values{1995}0.0000.0590.0000.0000.0560.0000.0590.000
process[Recruitment].ycs_values{1996}0.000-0.0150.0020.002-0.0080.000-0.0150.002
process[Recruitment].ycs_values{1997}0.0000.018-0.001-0.0010.0150.0000.018-0.001
process[Recruitment].ycs_values{1998}0.0000.014-0.000-0.0000.0130.0000.014-0.000
process[Recruitment].ycs_values{1999}0.000-0.0890.0000.000-0.0890.000-0.0890.000
process[Recruitment].ycs_values{2000}0.000-0.0350.0000.000-0.0340.000-0.0350.000
process[Recruitment].ycs_values{2001}0.000-0.0240.0000.000-0.0240.000-0.0240.000
process[Recruitment].ycs_values{2002}0.0000.0800.0000.0000.0800.0000.0800.000
process[Recruitment].ycs_values{2003}0.0000.3290.0000.0000.3290.0000.3290.000
process[Recruitment].ycs_values{2004}0.0000.4530.0000.0000.4530.0000.4530.000
process[Recruitment].ycs_values{2005}0.0000.4080.0000.0000.4080.0000.4080.000
process[Recruitment].ycs_values{2006}0.0000.4400.0000.0000.4400.0000.4400.000
process[Recruitment].ycs_values{2007}0.000-1.4990.0000.000-1.4990.000-1.4990.000
process[Recruitment].ycs_values{2008}0.0000.0930.0000.0000.0930.0000.0930.000
process[Recruitment].ycs_values{2009}0.0000.0540.0000.0000.0540.0000.0540.000
process[Recruitment].ycs_values{2010}0.0000.0300.0000.0000.0300.0000.0300.000
process[Recruitment].ycs_values{2011}0.000-0.1110.0000.000-0.1110.000-0.1110.000
process[Recruitment].ycs_values{2012}0.000-0.0540.0000.000-0.0540.000-0.0540.000
process[Recruitment].ycs_values{2013}0.000-0.0620.0000.000-0.0620.000-0.0620.000
process[Recruitment].ycs_values{2014}0.0000.0010.0000.0000.0010.0000.0010.000

CASAL objective function component values
ComponentBase_ModelSensitivity_1
potCPUE-51.180-51.180
potSurveyAbd-3.254-3.254
logbookLF54.66054.660
potCommAF54.41054.410
potSurveyAF118.200118.200
recCatchLF68.04068.040
prior_on_initialization.B09.5599.559
prior_on_recruitment.YCS9.6589.658
prior_on_size_at_age.Linf_male0.0000.000
prior_on_size_at_age.cv2_male0.0000.000
prior_on_selectivity[potSsel].male0.0000.000
prior_on_selectivity[potSsel].female0.0000.000
prior_on_q_qPotCPUE-8.331-8.331
prior_on_q_potSurveyq-11.380-11.380
clp1F0.0000.000
clp2F0.0000.000
clp3F0.0000.000
YCS0.5660.566
Total241.000241.000
Casal2 objective function component values
rownamesbetadiff_casal_flags_onbetadiff_casal_flags_offbetadiff_casal_flags_on_low_tolcppad_casal_flags_oncppad_casal_flags_offadolc_casal_flags_onadolc_casal_flags_offadolc_casal_flags_on_low_tol
observation->potSurvey-2010                 -1.038-1.042-1.038-1.038-1.042-1.038-1.042-1.038
observation->potSurvey-2014                 -1.020-1.022-1.020-1.020-1.022-1.020-1.022-1.020
observation->potSurvey-2018                 -1.196-1.197-1.196-1.196-1.197-1.196-1.197-1.196
observation->potCPUE-1990                 -0.162-0.159-0.162-0.162-0.159-0.162-0.159-0.162
observation->potCPUE-1991                 1.0981.1021.0981.0981.1021.0981.1021.098
observation->potCPUE-1992                 -1.057-1.053-1.057-1.057-1.053-1.057-1.053-1.057
observation->potCPUE-1993                 -1.991-1.990-1.991-1.991-1.990-1.991-1.990-1.991
observation->potCPUE-1994                 -2.144-2.143-2.144-2.144-2.143-2.144-2.143-2.144
observation->potCPUE-1995                 -2.085-2.085-2.085-2.085-2.085-2.085-2.085-2.085
observation->potCPUE-1996                 -2.186-2.186-2.186-2.186-2.186-2.186-2.186-2.186
observation->potCPUE-1997                 -1.960-1.959-1.960-1.960-1.959-1.960-1.959-1.960
observation->potCPUE-1998                 -2.157-2.157-2.157-2.157-2.157-2.157-2.157-2.157
observation->potCPUE-1999                 -1.689-1.689-1.689-1.689-1.689-1.689-1.689-1.689
observation->potCPUE-2000                 -2.192-2.192-2.192-2.192-2.192-2.192-2.192-2.192
observation->potCPUE-2001                 -2.164-2.164-2.164-2.164-2.164-2.164-2.164-2.164
observation->potCPUE-2002                 -2.042-2.042-2.042-2.042-2.042-2.042-2.042-2.042
observation->potCPUE-2003                 -2.172-2.173-2.172-2.172-2.173-2.172-2.173-2.172
observation->potCPUE-2004                 -2.192-2.192-2.192-2.192-2.192-2.192-2.192-2.192
observation->potCPUE-2005                 -1.908-1.911-1.908-1.908-1.911-1.908-1.911-1.908
observation->potCPUE-2006                 -1.793-1.797-1.793-1.793-1.797-1.793-1.797-1.793
observation->potCPUE-2007                 -2.176-2.177-2.176-2.176-2.177-2.176-2.177-2.176
observation->potCPUE-2008                 -2.188-2.189-2.188-2.188-2.189-2.188-2.189-2.188
observation->potCPUE-2009                 -2.051-2.051-2.051-2.051-2.051-2.051-2.051-2.051
observation->potCPUE-2010                 -2.122-2.124-2.122-2.122-2.124-2.122-2.124-2.122
observation->potCPUE-2011                 -1.875-1.865-1.875-1.875-1.865-1.875-1.865-1.875
observation->potCPUE-2012                 -1.767-1.756-1.767-1.767-1.756-1.767-1.756-1.767
observation->potCPUE-2013                 -1.775-1.769-1.775-1.775-1.769-1.775-1.769-1.775
observation->potCPUE-2014                 -1.038-1.034-1.038-1.038-1.034-1.038-1.034-1.038
observation->potCPUE-2015                 -1.851-1.851-1.851-1.851-1.851-1.851-1.851-1.851
observation->potCPUE-2016                 -2.092-2.092-2.092-2.092-2.092-2.092-2.092-2.092
observation->potCPUE-2017                 -2.104-2.105-2.104-2.104-2.105-2.104-2.105-2.104
observation->potCPUE-2018                 -1.352-1.347-1.352-1.352-1.347-1.352-1.347-1.352
observation->AFpotSurvey-2010                 39.65939.75439.65939.65939.75439.65939.75439.659
observation->AFpotSurvey-2014                 39.37139.38339.37139.37139.38339.37139.38339.371
observation->AFpotSurvey-2018                 39.16539.13639.16539.16539.13639.16539.13639.165
observation->AFpotFishing-2018                 47.10347.14647.10347.10347.14647.10347.14647.103
observation->AFpotFishing-2019                 7.8887.8877.8887.8887.8877.8887.8877.888
observation->lgobookLF                 54.89454.86254.89454.89454.86254.89454.86254.894
observation->rec.catch.LF              67.72768.61167.72767.72768.61167.72768.61167.727
prior->B0->process[Recruitment].b0              9.5619.5629.5619.5619.5629.5619.5629.561
prior->male_Linf->age_length[asMm0].linf            0.0000.0000.0000.0000.0000.0000.0000.000
prior->male_cv2->age_length[asMm0].cv_last         0.0000.0000.0000.0000.0000.0000.0000.000
prior->potSurvey_mA50->selectivity[potSurveySel_male].a50             0.0000.0000.0000.0000.0000.0000.0000.000
prior->potSurvey_mAto95->selectivity[potSurveySel_male].ato95           0.0000.0000.0000.0000.0000.0000.0000.000
prior->potSurvey_fA50->selectivity[potSurveySel_female].a50             0.0000.0000.0000.0000.0000.0000.0000.000
prior->potSurvey_fAto95->selectivity[potSurveySel_female].ato95           0.0000.0000.0000.0000.0000.0000.0000.000
prior->YCS->process[Recruitment].ycs_values{1980}1.2731.2751.2731.2731.2751.2731.2751.273
prior->YCS->process[Recruitment].ycs_values{1981}1.3151.3171.3151.3151.3171.3151.3171.315
prior->YCS->process[Recruitment].ycs_values{1982}1.3011.3021.3011.3011.3021.3011.3021.301
prior->YCS->process[Recruitment].ycs_values{1983}1.1841.1841.1841.1841.1841.1841.1841.184
prior->YCS->process[Recruitment].ycs_values{1984}0.9180.9170.9180.9180.9170.9180.9170.918
prior->YCS->process[Recruitment].ycs_values{1985}0.4770.4760.4770.4770.4760.4770.4760.477
prior->YCS->process[Recruitment].ycs_values{1986}-0.027-0.027-0.027-0.027-0.027-0.027-0.027-0.027
prior->YCS->process[Recruitment].ycs_values{1987}-0.304-0.304-0.304-0.304-0.304-0.304-0.304-0.304
prior->YCS->process[Recruitment].ycs_values{1988}-0.099-0.098-0.099-0.099-0.098-0.099-0.098-0.099
prior->YCS->process[Recruitment].ycs_values{1989}0.6870.6900.6870.6870.6900.6870.6900.687
prior->YCS->process[Recruitment].ycs_values{1990}1.0271.0241.0271.0271.0241.0271.0241.027
prior->YCS->process[Recruitment].ycs_values{1991}0.0820.0820.0830.0830.0820.0820.0820.083
prior->YCS->process[Recruitment].ycs_values{1992}-0.165-0.165-0.165-0.165-0.165-0.165-0.165-0.165
prior->YCS->process[Recruitment].ycs_values{1993}-0.166-0.166-0.166-0.166-0.166-0.166-0.166-0.166
prior->YCS->process[Recruitment].ycs_values{1994}0.1400.1410.1400.1400.1410.1400.1410.140
prior->YCS->process[Recruitment].ycs_values{1995}1.0061.0041.0061.0061.0051.0061.0041.006
prior->YCS->process[Recruitment].ycs_values{1996}0.1210.1210.1210.1210.1210.1210.1210.121
prior->YCS->process[Recruitment].ycs_values{1997}-0.212-0.212-0.212-0.212-0.212-0.212-0.212-0.212
prior->YCS->process[Recruitment].ycs_values{1998}-0.199-0.199-0.199-0.199-0.199-0.199-0.199-0.199
prior->YCS->process[Recruitment].ycs_values{1999}-0.307-0.307-0.307-0.307-0.307-0.307-0.307-0.307
prior->YCS->process[Recruitment].ycs_values{2000}-0.213-0.213-0.213-0.213-0.213-0.213-0.213-0.213
prior->YCS->process[Recruitment].ycs_values{2001}-0.262-0.261-0.262-0.262-0.261-0.262-0.261-0.262
prior->YCS->process[Recruitment].ycs_values{2002}-0.161-0.160-0.161-0.161-0.160-0.161-0.160-0.161
prior->YCS->process[Recruitment].ycs_values{2003}0.1870.1810.1870.1870.1810.1870.1810.187
prior->YCS->process[Recruitment].ycs_values{2004}-0.074-0.079-0.074-0.074-0.079-0.074-0.079-0.074
prior->YCS->process[Recruitment].ycs_values{2005}0.0920.0980.0920.0920.0980.0920.0980.092
prior->YCS->process[Recruitment].ycs_values{2006}-0.094-0.099-0.094-0.094-0.099-0.094-0.099-0.094
prior->YCS->process[Recruitment].ycs_values{2007}-0.285-0.279-0.285-0.285-0.279-0.285-0.279-0.285
prior->YCS->process[Recruitment].ycs_values{2008}-0.179-0.179-0.179-0.179-0.179-0.179-0.179-0.179
prior->YCS->process[Recruitment].ycs_values{2009}-0.058-0.058-0.058-0.058-0.058-0.058-0.058-0.058
prior->YCS->process[Recruitment].ycs_values{2010}0.5210.5210.5210.5210.5210.5210.5210.521
prior->YCS->process[Recruitment].ycs_values{2011}-0.294-0.294-0.294-0.294-0.294-0.294-0.294-0.294
prior->YCS->process[Recruitment].ycs_values{2012}-0.268-0.268-0.268-0.268-0.268-0.268-0.268-0.268
prior->YCS->process[Recruitment].ycs_values{2013}-0.286-0.286-0.286-0.286-0.286-0.286-0.286-0.286
prior->YCS->process[Recruitment].ycs_values{2014}-0.096-0.096-0.096-0.096-0.096-0.096-0.096-0.096
additional_prior->PrSurveyQ                 -11.367-11.368-11.367-11.367-11.368-11.367-11.368-11.367
additional_prior->PrCPUEQ                 -8.333-8.334-8.333-8.333-8.334-8.333-8.334-8.333
additional_prior->YCS_average_1                 0.5590.5590.5590.5590.5590.5590.5590.559
total_negloglike                 238.370239.365238.370238.370239.365238.370239.365238.370

## [1] "CASAL base model convergence information"
## [1] "Minimiser convergence threshold: 2e-005"                                                                 
## [2] "fmm: have converged: t = 1.53773e-005 f = 240.986"                                                       
## [3] "Successful convergence in optimise"                                                                      
## [4] "Minimiser achieved convergence after 91 quasi-Newton iterations using 101 objective function evaluations"
## [1] ""
## [1] "CASAL sensitivity 1 model convergence information"
## [1] "Minimiser convergence threshold: 1e-007"                                                                                 
## [2] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [3] "Successful convergence in optimise"                                                                                      
## [4] "Minimiser achieved convergence after 99 quasi-Newton iterations using 119 objective function evaluations"                
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model convergence information"
## [1] "fmm: have converged: t = 4.14863e-06 f = 238.37"
## [2] "Successful convergence in optimise"             
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model convergence information"
## [1] "fmm: have converged: t = 2.71219e-06 f = 239.365"
## [2] "Successful convergence in optimise"              
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model convergence information"
## [1] "fmm: step size too small. Indicates successful convergence (though this is not the textbook ideal convergence situation)"
## [2] "Successful convergence in optimise"                                                                                      
## [1] ""
## [1] "Casal2 cppad_casal_flags_on model convergence information"
##  [1] "Number of nonzeros in equality constraint Jacobian...:        0"
##  [2] "Number of nonzeros in inequality constraint Jacobian.:       42"
##  [3] "Number of nonzeros in Lagrangian Hessian.............:      903"
##  [4] "Number of Iterations....: 20"                                   
##  [5] "Number of objective function evaluations             = 110"     
##  [6] "Number of objective gradient evaluations             = 21"      
##  [7] "Number of equality constraint evaluations            = 0"       
##  [8] "Number of inequality constraint evaluations          = 110"     
##  [9] "Number of equality constraint Jacobian evaluations   = 0"       
## [10] "Number of inequality constraint Jacobian evaluations = 21"      
## [11] "Number of Lagrangian Hessian evaluations             = 20"      
## [12] "EXIT: Maximum Number of Iterations Exceeded."                   
## [13] "Number of nonzeros in equality constraint Jacobian...:        0"
## [14] "Number of nonzeros in inequality constraint Jacobian.:       42"
## [15] "Number of nonzeros in Lagrangian Hessian.............:      903"
## [16] "Number of Iterations....: 20"                                   
## [17] "Number of objective function evaluations             = 146"     
## [18] "Number of objective gradient evaluations             = 21"      
## [19] "Number of equality constraint evaluations            = 0"       
## [20] "Number of inequality constraint evaluations          = 146"     
## [21] "Number of equality constraint Jacobian evaluations   = 0"       
## [22] "Number of inequality constraint Jacobian evaluations = 21"      
## [23] "Number of Lagrangian Hessian evaluations             = 20"      
## [24] "EXIT: Maximum Number of Iterations Exceeded."                   
## [1] ""
## [1] "Casal2 cppad_casal_flags_off model convergence information"
##  [1] "Number of nonzeros in equality constraint Jacobian...:        0"
##  [2] "Number of nonzeros in inequality constraint Jacobian.:       42"
##  [3] "Number of nonzeros in Lagrangian Hessian.............:      903"
##  [4] "Number of Iterations....: 20"                                   
##  [5] "Number of objective function evaluations             = 115"     
##  [6] "Number of objective gradient evaluations             = 21"      
##  [7] "Number of equality constraint evaluations            = 0"       
##  [8] "Number of inequality constraint evaluations          = 115"     
##  [9] "Number of equality constraint Jacobian evaluations   = 0"       
## [10] "Number of inequality constraint Jacobian evaluations = 21"      
## [11] "Number of Lagrangian Hessian evaluations             = 20"      
## [12] "EXIT: Maximum Number of Iterations Exceeded."                   
## [13] "Number of nonzeros in equality constraint Jacobian...:        0"
## [14] "Number of nonzeros in inequality constraint Jacobian.:       42"
## [15] "Number of nonzeros in Lagrangian Hessian.............:      903"
## [16] "Number of Iterations....: 20"                                   
## [17] "Number of objective function evaluations             = 154"     
## [18] "Number of objective gradient evaluations             = 21"      
## [19] "Number of equality constraint evaluations            = 0"       
## [20] "Number of inequality constraint evaluations          = 154"     
## [21] "Number of equality constraint Jacobian evaluations   = 0"       
## [22] "Number of inequality constraint Jacobian evaluations = 21"      
## [23] "Number of Lagrangian Hessian evaluations             = 20"      
## [24] "EXIT: Maximum Number of Iterations Exceeded."                   
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model convergence information"
## character(0)
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model convergence information"
## character(0)
## [1] "CASAL model warnings"
## NULL
## [1] ""
## [1] "CASAL sensitivity 1 model warnings"
## NULL
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 betadiff_casal_flags_off model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 betadiff_casal_flags_on_low_tol model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 cppad_casal_flags_on model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 cppad_casal_flags_off model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 adolc_casal_flags_on model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 adolc_casal_flags_off model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"
## 
## [1] ""
## [1] "Casal2 adolc_casal_flags_on_low_tol model warnings"
## $warnings_found
## [1] 1
## 
## $warning_0
## [1] "Estimates were removed because of matching lower and upper bounds. Originally had 127 estimates, now have 42"
## 
## $type
## [1] "warnings"

1.5 Matching of outputs

Time series comparisons with CASAL base model results

## [1] "Catch time series base model comparison for run betadiff_casal_flags_on"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_off"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run betadiff_casal_flags_on_low_tol"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run cppad_casal_flags_on"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run cppad_casal_flags_off"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_off"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""
## [1] "Catch time series base model comparison for run adolc_casal_flags_on_low_tol"
## [1] "Removals / Actual catches for line match: no"
## [1] "Discarded catches for line match: no"
## [1] "Removals / Actual catches for pot match: no"
## [1] "Discarded catches for pot match: no"
## [1] "Removals / Actual catches for rec match: no"
## [1] "Discarded catches for rec match: no"
## [1] ""

Derived quantities

SB0, SBcurrent, MSY, F_MSY, others...

1.6 Plots

Comparison plots

## [1] "CASAL base parameter correlation range (excluding 1.0): -0.6876 0.9999"

## [1] "CASAL sensitivity 1 parameter correlation range (excluding 1.0): -0.6894 0.9999"

## [1] "Casal2 betadiff_casal_flags_on parameter correlation range (excluding 1.0): -0.704685923556288 0.999820714207084"

## [1] "Casal2 betadiff_casal_flags_off parameter correlation range (excluding 1.0): -0.697745743037247 0.999864061627529"

## [1] "Casal2 betadiff_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.719268887862089 0.999845946100166"

## [1] "Casal2 cppad_casal_flags_on parameter correlation range (excluding 1.0): -0.708287635061791 0.678193985161796"

## [1] "Casal2 cppad_casal_flags_off parameter correlation range (excluding 1.0): -0.705974195067 0.685973588521794"

## [1] "Casal2 adolc_casal_flags_on parameter correlation range (excluding 1.0): -0.704683609734695 0.99981954195823"

## [1] "Casal2 adolc_casal_flags_off parameter correlation range (excluding 1.0): -0.697745743037247 0.999863100354264"

## [1] "Casal2 adolc_casal_flags_on_low_tol parameter correlation range (excluding 1.0): -0.719263523071944 0.999846316745602"